Large orbital magnetic moments in carbon nanotubes generated by resonant transport 
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The nonequilibrium Green's function method is used to study the ballistic transport in metallic 
carbon nanotubes when a current is injected from the electrodes with finite bias voltages. We reveal, 
both analytically and numerically, that large loop currents circulating around the tube are induced, 
which come from a quantum-mechanical interference and are much larger than the current along 
the tube axis when the injected electron is resonant with a time-reversed pair of degenerate states, 
which are, in fact, inherent in the zigzag and chiral nanotubes. The loop current produces large 
orbital magnetic moments, making the nanotube a molecular solenoid. 

PACS numbers: 75.75.+a, 72.20.My, 73.23.Ad, 73.63.-b 
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INTRODUCTION 

Carbon nanotubes (CNTs) have remarkable elec- 
tronic properties related with their unique geometrical 
structure.— ^ For example, the way in which the graphene 
sheet is wound into a cylinder determines whether CNTs 
are metals or semiconductors Recent advances have 
made the measurements of electrical properties of indi- 
vidual single-wall CNTs possible.— 1 ^ Most recently, or- 
bital magnetic moments of a single-wall CNT have been 
detected, where the presence of moment is deduced from 
the shift of energy levels in external magnetic fields.—^ 
Since the magnitude of magnetic moments estimated in 
the experiment is about ten times the Bohr magneton, 
the CNT as a molecular solenoid has attracted attention. 

Semiclassically, the large orbital magnetic moments 
can be understood as an effect of chiral currents, for 
which several semiclassical calculations have been done 
based on Boltzmann's equation ! 10 : 11 However, the ap- 
proach, being applicable when the transport is diffu- 
sive and lacks coherence, cannot treat purely quantum 
mechanical effects. Quantum mechanical treatments 
of magnetic properties of CNTs have previously been 
given ) 12 ' 13 i 14 ' 15 but those studies have concentrated on 
the energy shift and moment of an isolated CNT in ex- 
ternal magnetic fields in equilibrium, so the currents and 
the effects of electrodes have not been studied. Since 
ballistic transport can be strongly affected by electrodes, 
we need to take account of the influence of electrodes in 
CNTs, which is exactly our motivation here. In the pro- 
cess, the quantum loop current, a purely quantum me- 
chanical effect, has turned out to be indeed large, which 
we study by making use of the nonequilibrium Green's 
function method^ 

A loop current has been studied in the context of small 
molecules placed between the scanning tunneling micro- 
scope (STM) tip and the substrate. Specifically, Nakan- 
ishi and TsukadaA£ and Tsukada et a/.— have shown that 
a loop current can be dramatically amplified when the 
energy of the injected electron is resonant with degener- 
ate eigenstates of an isolated molecule. This mechanism 
should also be applicable to CNTs, for which a special 
interest is the effect of inherent degeneracies associated 



with wave functions traveling clockwise and anticlockwise 
around the tube. We develop a perturbation expansion 
of the nonequilibrium Green's function in the weak cou- 
pling between the conductor and electrodes. This serves 
to identify how the loop current paths are determined in 
terms of phase variation of wave functions. The peculiar- 
ity of the CNTs appears as their band structure that can 
satisfy the condition for large orbital magnetic moments. 

We then numerically calculate the current distribution 
and the magnitude of the orbital magnetic moments as a 
function of the bias voltage, which confirms the analytic 
discussions. The order of magnitude of the magnetic mo- 
ment obtained in the numerical calculation is consistent 
with the value estimated in the experiment.— Through- 
out this Brief Report, we consider metallic CNTs in the 
absence of external magnetic fields. CNTs are character- 
ized by the chiral index (n\, n-i) in standard literatures,-^ 
where CNTs are metals when ri\ — ri2 = (mod3) or 
semiconductors otherwise^^ Here we concentrate on the 
one-body problem. 



FORMALISM 



We employ the tight-binding model whose Hamiltonian 
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Here, Hqnt is the tight-binding Hamiltonian for the 
single- wall CNT, -ffcNT-oicctrodc describes the connection 
between the CNT and the electrodes, ^electrode describes 
the electrodes, s,s',d,d' are the sites connecting the CNT 
and the electrodes and t,t',t"(> 0) are the respective 
transfer integrals [see Fig. HJa)]. t, the transfer integral 
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FIG. 1: (Color online) (a) A nanotube attached to ID elec- 
trodes is schematically shown. The arrows represent the quan- 
tum loop current, (b) A typical asymmetric peak of the cur- 
rent density Jy vs energy. 



with E in = i(Ef - Y, A ) at zero temperature These 
currents generate an orbital magnetic moment M = 
±fdrrx j(r) = \ I^n x r r 

Let us consider the case where the incident electron 
is resonant to doubly degenerate states 4>i ; ^2 with an 
eigenenergy e a . In the degenerate perturbation theory, 
the zeroth-order <f) R (= tfi A ) and 4> R {— 4*2)1 having 
E\,2 = £q + Ef 2 with Ef 2 the eigenvalues of 2 x 2 T, R , 
can be written as {<p\ ± (f>2)/V2, respectively, where we 
fix <pi 2 as the time-reversed pair with <pi = cfy^ = 4> a and 
choose the overall phase of <f> a so that [</> Q ( s )] 2 + [4>a{d)] 2 
be real. The Green's function becomes 



G R (i,j) 
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Im[(/) Q («)]Im[0 a (j)] 
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of CNT, is 0(1 eV)±, while we have taken t" for the elec- 
trodes to be large enough so that their density of states 
is nearly constant (wide-band limit) and t' , the hopping 
between CNT and electrodes, to be <C t, t" . Here, we 
follow other theoretical literature in assuming electrodes 
to be one dimensional for simplicity. 

For one-dimensional (ID) electrodes, the retarded and 
advanced Green's functions G R ' A (E) = (E - H ± ie)~ x 
[where + (— ) is assigned to R (A)] can be cast into a 
form— 
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4 c p (p = s,d), 



where A = (t') 2 /(tt") characterizes the sample-electrode 
coupling, and k ln the Bloch's wave number for the inci- 
dent electron having an energy E = —2t" cos(fci n ao) with 
ao the lattice constant of the conductor and electrodes. 
This reduces the problem to the sample alone, with the 
effects of electrodes contained in the self-energy T, R ' . 

We can perform perturbation expansions in terms of 
the weak coupling A. Since nonperturbative effects of- 
ten arise in resonance phenomena, here we adopt the 
eigenfunction expansion method^ that can take non- 
perturbative corrections into account. In this formal- 
ism, we express the Green's function as G R (ri,r2) = 
Ea^a(n)<^>2)/(£ - E a ), where (j> R > A {r) are the 
eigenfunctions of the non-Hermitian operators -ffcNT + 
T, R ' A that form a biorthonormal set. The current Zy flow- 
ing from site i to j is given by ~ f * dEJij (E) , where 
ix s and Hd are the chemical potentials of the source and 
drain electrodes, respectively. The current density Jy 
per unit energy is calculated as 



J ij (E) = ^Jm[H ij G] i (E% 



(4) 



where G n is the nonequilibrium Green's function de- 
fined by the usual Green's functions as G n = G R T, ln G A 



Here, we have only retained the term related to the res- 
onant states cv, which is valid as long as the splitting, 
~ At, of the degenerate levels due to the self-energy is 
smaller than the interval across the adjacent energy level 
[~ 1 meV/L(^m) for CNTs£], i.e., A < 1(T 3 /£(a™) for 
CNTs with t ~ 1 eV. 

From Eqs. (|4|) and ([5]), we end up with 
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This expression can also be derived by using the results 
of Nakanishi and TsukadaJi We can see that the E de- 
pendence of the current density is mainly determined 
by the factor (E - e a )/{\E -s a - Z R \ 2 \E - s a - T, R \ 2 ). 
This function exhibits an asymmetric peak [with width 
~ A and height ~ 1/A, see Fig. QXb)], which is why the 
integration across the peak results in a finite net cur- 
rent. Sasada and Hatano^ has interpreted the asym- 
metric resonance shape for the transmission coefficient 
as Fano effect, i.e., an interference between continuous 
states in electrodes and discrete states in a conductor. 
So the present case may be called a Fano effect extended 
to loop currents. 



ANALYTIC EXPRESSION FOR THE CURRENT 
AND MAGNETIC MOMENT 

When we integrate J z " (E) to obtain /PS , we can replace 
the range fid < E < fi s with — 00 < E < 00, since 
J!j* is negligible outside the peak. By taking care of E 
dependences in h ln and Ef 2 as well, we obtain 
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FIG. 2: (Color online) Band structure against k± for typical 
(a) zigzag or (b) armchair CNTs. The vertical dashed lines 
indicate discrete k± points. The horizontal line represents the 
energy of the incident electron. 




Since T, R ~ A, the leading term of is ~ A . This is to 
be contrasted with a nondegenerate resonant state a, for 
which the corresponding expression involves the factor 
Im[(/> Q (z)0* (J)], but this must vanish in the absence of 
external magnetic fields, so the current is small [0(A)] . 
For the degenerate case, /?■ oc lm[0 a (i)^>* (J)] indicates 
that currents flow along the direction in which the phase 
of <f> a varies. So this implies that a cylindrical conductor, 
with a standing wave along the axis and a propagating 
wave around the tube, should have currents circulating 
around the tube that is much larger than the component 
flowing along the axis. 

From Eq. ©, the magnetic moment along the tube 
axis generated by the current becomes 

(M , 2e£ lm{[r a (s)Md)?} 
[iviah - h t // £ «[|^( s )|2 + |^ (d) |2 ]2 

x £ln# a (i)C(.7)K^ x rj) z +0(X). (7) 

(ij) 

The total magnetic moment is the sum over the reso- 
nant states which are relevant to transport. While M a 
is proportional to the resonant energy level e Q , this is 
only valid for ID electrodes with Eq. so the specific 
form should depend on the detail of the electrodes. In 
addition, we note that, while one might expect that edge 
states, which are known to exist in a finite carbon nan- 
otube with a flat dispersion on the Fermi energy^ may 
have an important effect on the magnetic moment, the 
contribution from edge states is negligible for ID elec- 
trodes since M a shown in Eq. (J7|). 

Equation ([7]) contains a factor lm{[0* (s)(f> a (d)] 2 }, 
which gives M a cx sin(2fc"i?y>) if we assume, for a heuris- 
tic purpose, a simple plane- wave form around the circum- 
ference 4> a {s) = <j)* a (d) — \(j> a (s)\exp(ik"Rip/2), where 
kj_ is the wave number around the tube of the state a, R 
is the radius of the tube, and <p is the angle subtended by 
the two electrodes, as shown in Fig. [IJa). Only when the 
electrodes are attached asymmetrically (<p ^ 0) does the 
loop current arise, which immediately resolves a puzzle 
one might have on the symmetry: how a specific sense 



FIG. 3: (Color online) (a) Typical doubly degenerate (in black 
and red) eigenfunctions with an eigenenergy of 0.419 522t in a 
zigzag CNT with (ni,n 2 ) = (18,0) and L — Ua /V3. o and* 
indicate the sign and the amplitude of the wave function. The 
dashed line indicates how the honeycomb lattice is wound into 
a tube, (b) A typical current distribution which is resonant 
to the state illustrated in (a). The size of the colored arrows 
indicates the magnitude of the current density, while the black 
arrows indicate the source and drain electrodes. 



of rotation of the circulating current can arise when the 
cylindrical conductor has no structural chirality (as in 
zigzag or armchair CNTs). 

The above expression for M contains another fac- 
tor, Im[0 Q (i)0* (j)] tx sin{fc"[#(i) — 6(j)]}, where we set 
4>a(r) = |</> Q ,(r)|e ifej - e ( T '' ) with 9 the coordinate along the 
circumference. This factor expresses that loop currents 
are generated due to interference between the doubly de- 
generate states. Since these resonant states (encircled in 
Fig. [2| that participate in electronic transport lie near 
Ep — 0, the range of the value of k± is determined by the 
band structure of the conductor in general. In particular, 
the metallic CNTs with {n\ — 71-2) = (mod 3) have 2 - (see 
Fig. HI) 

JO, (ni —ri2)jn = (mod 3) : armchair CNT 
- 1 = \ iff , (m - n 2 )/n ^ (mod 3) : zigzag CNT, 

where n is the greatest common divisor of (rii, 712) and a 
the length of the translational vector along the circum- 
ference. Thus, armchair CNTs with k± = should have 
nearly zero magnetic moments, which is contrasted with 
zigzag CNTs for which magnetic moments are large. Chi- 
ral CNTs can have either {nx — n-z)/n = or ^ (mod 3), 
so even chiral structures, which may naively seem to 
guarantee a molecular solenoid, can have large moments 
only when the dispersion is correct. 

So we conclude that large loop currents arise when 
we have (i) time-reversed pair of degenerate states, (ii) 
asymmetrically attached source and drain electrodes, and 
(iii) k± + 0. 
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FIG. 4: (Color online) Calculated total current flowing (a) 
along the axis and (b) around the tube versus the bias voltage 
in a zigzag CNT, for which we have chosen smaller (m, ri2) = 
(3, 0) for clarity with L = 20ao /v3 with the same asymmetry 
(Rip) in the electrodes as in Fig. 0R). Note the different 
vertical scales. 



M a ~ (e/h)e a RaQ. Around the Fermi energy, a CNT 
has a cone-like dispersion with e a = [1 mcV/L (/mi)] x 
a (a = 1,2,...). Summing over the energy levels be- 
tween fj, s = Vbias and fid = 0, we obtain M — J2 a M a ~ 
10 _3 //_bxL (/xm) x [Vbias (meV)] 2 x(i?/a ). This amounts 
to M ~ 10[ib for L ~ l ;/xm, Vbias ~ 30 meV, and 
R ~ 10ao. We also consider the magnetic field B ln< \ 
generated by the circulating currents. We must sat- 
isfy M ■ B ind = M(fi /2TrR)(M/nR 2 ) < At for the 
energy shift due to this to be negligible. If we com- 
bine this with the previous condition, we have to have 
[L (fimtfiR/ao)- 1 [V hias (eV)] 4 « A « 10- 3 [L HI" 1 
for t ~ 1 eV. For Vbi as and R assumed above, the above 
set of inequalities is satisfied for L < 1 /im. 



DISCUSSION 



NUMERICAL RESULTS FOR CARBON 
NANOTUBES 

We now move on to the numerical calculation to con- 
firm the analytical discussions above. We have calculated 
the current and the magnetic moment for various types 
of CNT directly from Eq. © for t'/t = 0.2 and t" /t = 2 
with A = 0.02. Figure[3Ja) shows typical doubly degener- 
ate eigenstates in a zigzag CNT, while the current density 
resonant to the state is depicted in Fig. [Hb). We can 
immediately sec that the current circulating around the 
tube is much larger than the one flowing along the axis, 
which endorses the discussion above. 

Figure [4fa) shows the total current flowing along the 
axis of the CNT as the function of the bias voltage. We 
can see that the current along the axis is nearly dis- 
cretized, which is due to the resonant tunneling from a 
ID electrode through the discrete electron levels existing 
at regular intervals in a CNT. FigurefJJb) shows the total 
current circulating around the tube versus bias voltage, 
which is indeed larger than the current along the axis 
with a ratio ~ A -1 S> 1. The behavior consists of steps, 
but this time the increment of the current is not constant. 
This is attributed to the fact that the contribution to the 
current (or the orbital magnetic moment) from a pair of 
resonant eigenstates a differs from pair to pair [oc e a , Eq. 
©, for ID electrodes]. 

EVALUATION OF THE MAGNETIC MOMENT 

Let us finally estimate the order of magnitude of the 
orbital magnetic moment in CNTs. From Eq. ijTj). 



Here, we have concentrated on the one-body problem. 
The effects of the electron-electron interaction appear in, 
e.g., Coulomb blockade^ and Tomonaga-Luttinger liq- 
uid stated So we need to take the interaction into ac- 
count to describe the real systems more accurately, which 
is a future problem. Another point is that, although the 
self-induced magnetic field is small as mentioned above, 
some CNTs have large magnetic susceptibility, * 2 so its 
effect may be an interesting problem. Also assumed is 
that each electrode touches only one atom of the CNTs. 
We have checked that the loop current is not significantly 
affected even when each electrode touches more than one 
atoms in the CNT. 

Although orbital magnetic moments of CNTs come to 
be detected indirectly from the shift of energy levels in 
external magnetic fields, some direct experimental obser- 
vation is desirable. We may be able to use the STM 22 
to probe the loop current or attach electrodes to CNT 
ropes or "forest"^ to detect the overall magnetization. 
For the latter, some statistical average may have to be 
involved. 
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